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In order to obtain stable and accurate general relativistic simulations, re-formulations of the Ein- 
stein equations are necessary. In a series of our works, we have proposed using eigenvalue analysis of 
constraint propagation equations for evaluating violation behavior of constraints. In this article, we 
classify asymptotical behaviors of constraint-violation into three types (asymptotically constrained, 
asymptotically bounded, and diverge), and give their necessary and sufficient conditions. We find 
that degeneracy of eigenvalues sometimes leads constraint evolution to diverge (even if its real-part 
is not positive), and conclude that it is quite useful to check the diagonalizability of constraint 
propagation matrices. The discussion is general and can be applied to any numerical treatments of 
constrained dynamics. 



I. INTRODUCTION 

So-called numerical relativity (computational simula- 
tions in general relativity) is a promising research field 
having implications for ongoing astrophysical observa- 
tions such as gravitational wave astronomy [1]. Many 
simulations of binary compact objects have revealed that 
mathematically equivalent sets of evolution equations 
show different numerical stability in the free-evolution 
scheme. 

There are many approaches to re-formulate the Ein- 
stein equations for obtaining a longterm stable and ac- 
curate numerical evolution (e.g. see references in [2]). In 
a series of our works, we have proposed the construction 
of a system that has its constraint surface as an attrac- 
tor. By applying eigenvalue analysis of constraint prop- 
agation equations, we showed that there is a constraint- 
violating mode in the standard Arnowitt-Deser-Misncr 
(ADM) evolution system [3, 4] when it is applied to a sin- 
gle non-rotating black- hole space-time [6]. We also found 
that such a constraint-violating mode can be compen- 
sated for if we adjust the evolution equations with a par- 
ticular modification using constraint terms like the one 
proposed by Detweiler [5] . 

Our predictions are borne out in simple numerical ex- 
periments using the Maxwell, Ashtekar, and ADM sys- 
tems [6, 7, 8, 9]. There are also several numerical exper- 
iments to confirm our predictions are effective[10, 11]. 
However we have not yet obtained definite guidelines for 
specifying the above adjusting terms and their multipli- 
ers. 

In this article, we show the essential steps in analyz- 
ing constraint amplification factors (defined in §11 B). In 
§111, we show that degeneracy of eigenvalues sometimes 
leads constraint evolution to diverge. This observation 
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suggests the importance of checking the diagonalizabil- 
ity of characteristic matrices, and gives further insights 
for constructing an asymptotically constrained system. 

II. A GUIDELINE TO OBTAIN A ROBUST 
EVOLUTION SYSTEM 

A. Idea of Adjusted system 

We begin by reviewing our proposal for an "adjusted 
system" . 

Suppose we have a dynamical system of variables 
u a {x l ,t), which has evolution equations, 

d t u a = f(u a ,d lU a ,---), (2.1) 

and the (first class) constraints, 

C a (u a ,d iU a ,---)^0. (2.2) 

Note that we do not require (2.1) to form a first-order 
hyperbolic form. We propose to investigate the evolution 
equation of C a (constraint propagation), 

d t C a =g{C a ,d % C a , •••), (2.3) 

for evaluating violation features of constraints. 

The character of constraint propagation, (2.3), will 
vary when we modify the original evolution equations. 
Suppose we modify (adjust) (2.1) using constraints 

d t u a = f(u a , d t u\ ■■■)+ Ftr.ftC", • • •), (2.4) 

then (2.3) will also be modified as 

d t C a = g{C a ,d t C a , ■■■) + G(C a , d l C a , ■ ■ •). (2.5) 

Therefore, finding a proper adjustment F(C a , • • •) is a 
quite important problem. 

Hyperbolicity analysis may be a way to evaluate con- 
straint propagation, (2.3) and (2.5) [12]. However, this 
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requires (2.3) to be a first-order system which is easy to 
be broken. (See e.g. Detweiler-type adjustment [5] in the 
ADM formulation [6] ) . Furthermore hyperbolicity analy- 
sis only concerns the principal part of the equation, that 
may fail to analyze the detail evaluation of evolution. 

Alternatively, we have proceeded an eigenvalue analy- 
sis of the whole RHS in (2.3) and (2.5) after a suitable 
homogenization, which may compensate for the above 
difficulties of hyperbolicity analysis. 

B. CP matrix and CAF 

We propose to transform the constraint propagation 
equation, (2.3) and (2.5), into Fourier modes, 

d t C a = g(C a ) = M a f3 CP, 
where C{x, t) a = j C(k, t) a cxp{ik ■ x)d 3 k, (2.6) 

then to analyze the eigenvalues, say A Q , of the coeffi- 
cient matrix, M a fj. We call A a and M a p the constraint 
amplification factors (CAFs) and constraint propagation 
matrix (CP matrix), respectively. 

So far we have proposed the following heuristic conjec- 
tures [6, 7, 8, 9]: 

(A) If the CAF has a negative real-part (the constraints 
are necessarily diminished), then we see more stable 
evolution than a system which has a positive CAF. 

(B) If the CAF has a non-zero imaginary-part (the con- 
straints are propagating away), then we see more 
stable evolution than a system which has a zero 
CAF. 

We observe that this eigenvalue analysis requires the 
fixing of a particular background space-time, since the 
CAFs depend on the dynamical variables, u a . 

C. Classification of Constraint propagations 

The CAFs indicate the evolution of constraint viola- 
tions (definitely its Fourier modes). It is natural to as- 
sume that a divergence of constraint norm is related to 
the numerical blow-ups. Therefore we classify the fun- 
damental evolution property of constraint propagation 
equation (2.6) as follows: 

(CI) Asymptotically constrained : Violation of con- 
straints decays (converges to zero). 

(C2) Asymptotically bounded : Violation of constraints 
is bounded at a certain value. 

(C3) Diverge : At least one constraint will diverge. 

Note that (CI) c (C2). We will derive the necessary 
and sufficient conditions for (CI) and (C2) in the next 
section. 



III. CONDITIONS FOR (CI) AND (C2) 

A. Preparation 

Hereafter, we consider a set of evolution equations, 

d t C l {t) = M l 3 C\ (3.1) 

where C l (i — 1, • • • , n) is a complex- valued vector, M % j 
is a n x n complex- valued matrix, and C l (t) is assumed 
to have finite- valued initial data C J (0). 

Without loss of generality, the CP matrix M can be as- 
sumed to be a Jordan normal form, since within complex- 
valued operations all the matrices can be converted to 
this form. Suppose that M has r different eigenvalues 
(Ai, • • • , A r ), where r < n. Let the multiplicity of Afe as 
n k , where Y^k=i n k — n. M can be expressed as 

• • °\ 

M = J 1 +...+J r : = , (3.2) 

\0 J r J 

where the cell size of J k is n k x rife. The Jordan matrix 
Jk is then expressed using a Jordan block J k e, 

Jk = Jki+ ■ ■ ■ +Jkm, (3.3) 

/A fe 1 0\ 
Jki := ■■. ■■. • (3.4) 

V O A fe 1 / 

Note that J kl is n M x n ke , Y^i=\ n kl = n k , m = n - 
rank(M — X k E) and max^nw) = v k . The minimum 
polynomial of M is written as 

fi M (t) - (t-Ai)" 1 ■■■(*- \ r ) Vr . (3.5) 

If Jfe is diagonal (i.e. n k = n — rank(M — \ k E)), then 
v k = 1 for that k. If M is diagonalizable (i.e. n k = 
n — rank(M — Afe-E) for Vfc), then v k = \ for all k. 
We then have the following statement. 

Proposition 1 The solution of 

d t C a - J k C a (3.6) 

can be expressed formally as 

Vk — l 

C a (t) = exp(Afei) J2 («f )te )- ( 3 - 7 ) 

i=o 

A proof is available by mathematical induction. Suppose 
that Jfei is v k x v k which is the maximal size J k . By direct 
calculation, we have that d t C a = J k \C a yields (3.7) with 
t-polynomial of degree (v k — 1). Then we see that (3.7) 
is satisfied in general. 

From this proposition, the highest power A^fe in 
t— polynomial in (3.7) is bounded by < Afe < v k — 1. 
The matrix J k in (3.6) can be directly extended to the 
full CP matrix, M, in (3.1). Therefore the highest power 
N in all constraints is bounded by 

< N < max (n k ) - 1. (3.8) 

l<fe<r 
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B. Asymptotically Constrained CP 

The following Propositions 2 and 3 give us the next 
theorem. 

Theorem 1 Asymptotically constrained evolution (vio- 
lation of constraints converges to zero) is obtained if and 
only if all the real parts of the CAFs are negative. 

Proposition 2 All the real part of CAFs are negative => 
Asymptotically constrained evolution. 

proof) We use the expression (3.7). If Ke(Afc) < for V/c, 
then Ci will converge to zero at t —> oo no matter what 
the t— polynomial terms are. D 

Proposition 3 Asymptotically constrained evolution => 
All the real parts of the CAFs are negative. 

proof) We show the contrapositive. Suppose there exists 
an eigenvalue Ai of which the real-part is non-negative. 
Then we get d t C\ — \\C\ of which the solution is C\ — 
Ci(0) exp(Aii). C\ does not converge to zero. □ 



C. Asymptotically Bounded CP 

The following Propositions 4 and 5 give us the next 
theorem. 

Theorem 2 Asymptotically bounded evolution (all the 
constraints are bounded at a certain value) is obtained 
if and only if all the real parts of CAFs are not positive 
and Jk is diagonal when SRe(Afc) = 0. 

Corollary Asymptotically bounded evolution is obtained 
if the real parts of CAFs are not positive and the CP 
matrix M a p is diagonalizable. 

Proposition 4 All the real parts of CAFs are not posi- 
tive and Jk is diagonal when Ke(Afe) = => Asymptoti- 
cally bounded evolution. 

proof) We use the expression (3.7). When !>Re(Afc) < 0, 
exp(Afci) x (i-polynomials) will converge to zero no matter 
what the t— polynomial terms are. When 9?e(Afc) = 0, we 
see Vk = 1 from the assumption of diagonality of So 
we see the t— polynomial terms are constant and exp(Afci) 
is bounded. □ 

Proposition 5 Asymptotically bounded evolution => All 
the real parts of the CAFs are not positive and Jk is di- 
agonal when 5R(Afe) = 0. 

proof) We show the contrapositive. If there exists an 
eigenvalue of which the real-part is positive, then con- 
straints will diverge no matter what the t— polynomial 
terms are. Therefore we try to show that constraints will 
diverge when all the real-parts of eigenvalues are non- 
positive, and there exists Xk such that 5J(A/j) = and its 
Jordan matrix Jk is not diagonal. 



Since Jordan matrix Jk is not diagonal, we see the 
power of t— polynomial Vk is greater than 1 in the ex- 
pression (3.7). Thus we have that (3.7) will diverge in 
t -> oo. n 



IV. CONCLUDING REMARKS 

Two theorems will give us a guideline to analyze a 
constraint-violating mode of the system. The result sup- 
ports our previous heuristic conjecture (A), but also sug- 
gests an ill-behaving case when CAFs are degenerated 
and its real-part is zero, when the associated Jordan ma- 
trix is not diagonal. This indicates the importance of 
checking the diagonalizability of constraint propagation 
matrix M. 

Along the line of our evaluation of constraint propa- 
gation equations (3.1), we propose a practical procedure 
for this classification in Figure 1 . We think that this dia- 
gram will provide systematic predictions for obtaining a 
robust evolution system in any constrained dynamics. 

The present classification is only on the fixed back- 
ground spacetime and only for t — > oo. It is still not 
clear at what value the constraints are bounded if a lim- 
iting value exists. Thus further modifications are un- 
derway. We are also applying the present classification 
scheme to various adjusted systems of the Einstein equa- 
tions (adjusted ADM, and further modified versions), to- 
gether with numerical experiments. We hope to report 
on them in the near future. 

The current constraint analysis only concentrates to 
the evolution equations and does not include the effect of 
the boundary treatments. Since the eigenvalues are eval- 
uated locally, it will be possible to include the effect of 
numerical boundary conditions if they are expressed ap- 



Ql: Is there a CAF which real part is positive? 
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Q2: Are all the real parts of CAFs negative? 
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Q5: Is the associated Jordan matrix diagonal? 
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FIG. 1: A flowchart to classify the fate of constraint propa- 
gation. 
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parently in a part of the evolution equations. This is also 
the one direction to proceed our future research. Mean- 
while, we would like to remark that one of our proposed 
adjustments in [9] contributes to enforce the computa- 
tional ability of the black-hole excision boundary treat- 
ment [11]. 

By extending the notion of "norm" or "compactness" 
of constraint violations, it might be interesting to define 
a new measure which monitors a "distance" between the 
constraint surface and an evolution sector in constraint 
dynamics. 
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